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Electric quadrupole and magnetic octupole of the A 
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Using a covariant spectator constituent quark model we predict an electric quadrupole moment 
Q/^+ = —0.043 efm^ and a magnetic octupole moment 0^+ — —0.0035 efm"^ for the A"*" excited 
state of the nucleon. 



Although it was the first nucleon resonance to be dis- 
covered, the properties of the A are almost completely 
unknown. Only the A"*""*" and A"*" magnetic moments 
have been measured, and these measurements have large 
error bars [H, 0, Si- Most of the information we have 
about the A comes from indirect information, such as 
the study of the 7iV — > A transition [1] . 

The dominant A elastic form factors are the electric 
charge Geo and magnetic dipole Gmi- The subleading 
form factors are the electric quadrupole {Ge2) and mag- 
netic octupole (Cms)- Those form factors measure the 
deviation of the charge and magnetic dipole distribution 
from a symmetric form [S]. At = the form fac- 



Ga/i(0). 



and the magnetic 



tors define the magnetic dipole 
electric quadrupole Qa = Ge2{^)-[^ 
octupole Oa = G'Af3(0)2]gx moments, where e is the 
electric charge and Ma the A mass. 

Until recently, there were essentially only theoretical 
predictions for /za (see Ref. @ for details) and Qa 0) B 
B [13, HHj dj [13, El- The exception was the pioneering 
work in lattice QCD [l^ , where all the form factors were 
estimated for low Q^, although the statistics for Ge2 and 
Gmz were very poor. 

Recent lattice QCD calculations of all four form fac- 
tors over a limited range have revived interest in the A 
moments, especially the interesting quadrupole and oc- 
tupole moments [H, [O] ■ These results are obtained only 
for unphysical pion masses in the range of 350-700 MeV 
so some extrapolation to the physical pion mass is re- 
quired [is . il9i] . Still, in the absence of direct experimen- 
tal information, lattice QCD provides the best reference 
for theoretical calculations. Stimulated by these new lat- 
tice results the covariant spectator quark model J3J and 
chiral Quark-Soliton model (xQSM) have been used 
to estimate the A form factors. Simultaneously, a lattice 
technique based on the background-field method [2l| has 
been used to estimate the /xa with great precision [23 |. 
The octupole moment Oa has also been evaluated by 
Buchmann Q using a deformed pion cloud model, and 
QCD sum rules (QCDSR) have been used to estimate 



both Qa and Oa [23|. 

The size of the moments Qa and Oa tells us if the A 
is deformed, and in which direction. The nucleon, as a 
spin 1/2 particle, can have no electric quadrupole mo- 
ment [24| [although the possibility remains, as pointed 
out by Buchmann and Henley [2^, that it might be a 
collective state with an intrinsic quadrupole moment]. 
While the measurement of the quadrupole form factors 
for the — + A transition give s some information about 
the deformation of the A [2y|, it is very important to 
obtain an independent estimate [13, [13 • Motivated by 
these considerations, the Nicosia-MIT and the Adelaide 
groups are presently working on an evaluation of Ga/s us- 
ing lattice QCD [TtI . [28| . Also Ledwig and collaborators 
are working in the same subject (13 using the xQSM. 

In this Letter we use the covariant spectator formalism 
28 1 to evaluate Qa and Oa. Following previous work 
3C , [111 , we describe the A as a quark-diquark system 
composed of a S-state with an admixture of two D states 



D3 



6* 



(1) 



where a is the mixture coefficient of the D3 state [L = 2, 
S — 3/2) and h the mixture coefficient of the Dl state 
(L — 2, S = 1/2). Each of the states are separately 
normalized, so that N = l/Vl + a? + 6^. The S, Dl and 
D3 wave functions are products of spin-isospin (and, for 
the D states, L = 2) operators and an appropriate scalar 
wave function -05, iI)di and which depends only the 
square of the momentum [P — fc)^ of the off-shell quark, 
where k is the four-momentum of the on-shell diquark 

In this model 0, [H, [s^, [Ml, [H, [H the A current can 
be written as 



(2) 



where P_ (-P+) is the initial (final) A momentum, and 
the sum is over all polarizations (A) of the diquark, and 
the covariant integral = J 



where is the 



2 



diquark energy. Additional terms proportional to o?N'^, 
b^N^ and abN^ can be neglected if a and b are small. 
The quark current jj in Eq. ^ includes a dependence 
on the quark u and d charges and anomalous magnetic 
moments k„ and Kd- See Refs. @, [13 for details. 

The current ([2]) can be written in a standard form in- 
volving four basic form factors, denoted F*, i = 1 — 4. 
The electric and magnetic moments are linear combina- 
tions of these [1, IH, HI, [3^ , and at — 0, to first order 
in the mixing coefficients a and b, they become 

Geo(O) = N^caIs 

Gmi(0) = N^/aIs 

Ge2{0) = 3{aN^)eAlD3 

Gm3(0) = fAN^[aI'o3 + 2bJ'oi], (3) 

where /a = ga + Maka/Mn, 

eA = 5(l + T'3), KA = + K-Ta), 

K+ = 2Hu - Kd, K- = |k„ -I- ^Kd, (4) 

with T3 = diag(3, 1, —1, —3), and 

Im = lim- / bik,q,P+)^D3iP+,k)^siP-,k) 
r Jk 

^'di = lim- / b{k,q,P+)i'Di{P+,k)ijs{P-,k), 

with T =_Q^/(4M|^) and b{k,q,P+) « Y2oik) as defined 
in Ref. [30]. The S-state wave function is normalized to 
unity (so that X5 = 1), and to first order in the mixing 
coefficients a and b, N"^ — > 1 so Geo{^) — ga, giving the 
correct charge. The multipole moments E2 and M3 are 
fixed by the factors I'jji and 1^3, and are zero if there are 
no D states. In particular Ge2{Q) is determined only by 
X^3, although G'7\/3(0) can depend on a delicate balance 
between 1^3, ^'^^ coefficients a and b. 

To illustrate how lattice data can be used to constrain 
models, we show results from two different parameteriza- 
tions for the A wave functions. The first one, denoted by 
Spectator 1 (Sp 1), is model 4 of Ref. [13] ■ That model 
fixed the pion cloud contribution (using a simple param- 
eterization) and adjusted the remaining valence contri- 
bution to fit the 7iV — > A data. The second parame- 
terization, from Ref. [3l| and denoted Spectator 2 (Sp 
2), uses the same functional form for the valence part of 
the D-state wave functions, but fits the valence part of 
the wave function directly to the lattice data [36[. Be- 
cause the pion mass used in these lattice calculations is 
large, the pion cloud effects are negligible at the lattice 
"point" and provide a better determination of the va- 
lence quark contribution at that point. After the fit is 
made, the results are extrapolated to the physical " point" 
by replacing the masses of the nucleon, A, and p meson 
(all parameters that enter into the functional form of the 
wave functions and currents) to their physical masses. 
We believe that model Sp 2 gives a more reliable param- 
eterization of the A wave function, but we compare it to 
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-3.63 


-1.79 





1.79 


Buchmann (imp) [121 


-2.49 


-1.25 





1.25 


Buchmann (exc) [1^ 


-9.28 


-4.64 





4.64 
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-0.0226 





0.0226 
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±0.0057 
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Spectator 1 


-3.87 


-1.93 





1.93 


Spectator 2 


-3.36 


-1.68 





1.68 



Lattice: 



Quenched Wilson -0.81±0.29 
Dynamical Wilson -0.87±0.67 

Hybrid -"^-OGIiI'.II 
TAam 1: Summary ot existing theoretical and lattice results 
for Ge2{0)- Lattice data from Ref. [17[- Quenched Wilson 
has ttItt = 411 MeV, dynamical Wilson has = 384 MeV, 
and hybrid has m,r ~ 353 MeV. 



Ga/3(0) A++ A+ AO A- 

GP [5| -11.68 -5.84 5.84 

QCDSR [23] -0.0925 -0.0462 0.0462 

error ±0.0234 ±0.0117 ±0.0117 

Spectator 1 -0.046 -0.023 0.00084 0.024 

Spectator 2 -3.46 -1.70 0.063 1.82 

TABLE IL Summary of existing theoretical results for 
Gm3{0). GP stands for general parameterization (of QCD). 



model Sp 1 to show the impact of using the lattice data 
to constrain the fit. In the first model (Sp 1) there is a 
mixture of 0.88% of D3 state and 4.36% of Dl state; the 
second model (Sp 2) has a mixture of 0.72% for both the 
D3 and Dl states. 

In this letter we restrict our discussion to the moments 
Qa and Oa, which are extracted from the values of the 
form factors Ge2 and Gms at — 0. A more complete 
study will be presented in a future work [s^] ■ Our results 
are true predictions; once the 7^'^ — *■ A reaction has been 
described no additional parameters are adjusted. The 
results for Ge2{0) are presented in tableUand for Gm3{0) 
m table [ID These are obtained from the integrals I'e^ — 
-7.00 and = 1.59 for Sp 1 and 1^3 = -6.65 and 
X^i = 0.24 for Sp 2. 

Before comparing our results with other models note 
that Sp 1 and Sp 2 give similar predictions for the 
quadrupole moment but very different predictions for the 
octupole moment. Clearly the octupole moment is more 
sensitive to the details of the model, and it is only the 
strong constraint imposed by the lattice data that allows 
us to predict that 6*^3(0) ~ -1.70. 

Both tables compare our results with predictions of 
other models. In table U we include the classic nonrela- 
tivistic quark model (NRQM) from Isgur et al. 7] , where 
the tensor color hyperfine interaction requires a mixture 
of D-state quarks with S-state quarks. This description 



considers only the valence degrees of freedom, and the 
contribution for the electric quadrupole moment is de- 
termined by both the mixture coefhcients and a confine- 
ment parameter [1, [l2| . In these models the contribution 
for the electric quadrupole can be estimated in impulse 
approximation [13, from 



Mmp) 



2 2 



(5) 



where is the neutron squared radius in fm^ . Using a re- 
cent value of = —0.116 fm^, we obtain Qa+ — —0.0464 
fm^, or (Oi^ —1.81 in close agreement with the val- 
ues from Ref. [iO| quoted in the table. (For a review of 
the earlier results, see Ref. Similar results are ob- 

tained by Buchmann et al. |12l | using a constituent quark 
model with a D-state admixture [1, [11] with a slightly 
different confinement parameterization and an impulse 
approximation to the one-body current. 

In the same work 12] , an estimate of the nonvalence 
contributions, based on a two-body exchange current rep- 
resentative of the nonvalence degrees of freedom, is ob- 
tained. These nonvalence contributions are the dominant 
ones, and assuming no D-state admixture, can be esti- 
mated from 



(6) 



Although developed in the constituent quark formalism 
this relation is parameter independent [l2| . The ex- 
pression ^ has also been derived in the large Nc limit 
[13|. Later, the expression ^ was improved using a 
general parameterization (GP) of QCD [3, [H, [i^], with 
the inclusion of higher order terms, and used to extract 
(0) = -7.02 ± 4.05 from the 7iV ^ A electric 
quadrupole data 0] . All of these results seem to suggest 
that the contribution of the pion cloud to the quadrupole 
moment could be quite large. On the other hand, calcu- 
lations based on xPT [ll|, and recent results derived in 
a xQSM [13] all of which include the pion cloud, suggest 
that the pion cloud effect might be smaller than estimates 
based on Eq. From this we conclude that model cal- 
culations of the size of the pion cloud contribution to the 
quadrupole moment are inconclusive. 

Finally, the tables show the lattice QCD simulations 
[13] based on three different approaches: a quenched cal- 
culation using a Wilson action with u and d quarks, a 
dynamical calculation using a Wilson action including u 
and d sea quarks, and a hybrid action which also includes 
strange sea quarks. The lattice data is however limited by 
the significant error bars that prevent an accurate extrap- 
olation to = (assuming a dipole or an exponential 
dependence on Q^) and by heavy pion masses (which 
require an extrapolation in m-^). Even so, the size of the 
hybrid calculation may be an indicator that the meson 
cloud contribution to Ge2 is not negligible, although not 
comparable with ([6]). Quark models can be important 
for extrapolating the lattice data to = and to the 
physical pion mass. In any case, the predictions of our 
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FIG. 1: Ge2 and Gm3 form factors for A"*". The Ge2 lattice 
data is from Ref. [17| and Gmz lattice data is from Ref. [161. The 
lattice points for = are result of an extrapolation [ITII . 



model should be compared to other calculations of the 
valence quark contributions to these moments. 

The dependence of the A+ form factors Ge2 and 
Gm3 are shown in figure [TJ Our results are completely 
consistent with the dependence of the lattice calcula- 
tions [H, [l3]- Future lattice QCD simulations would be 
important for a more precise constraint on . 

In conclusion, using our best model (Sp 2) we predict 



Q 



A + 



-0.043 efm^ 



-0.0035 efm^ 



(7) 



This estimate for lies between the negligible predic- 
tions of QCD sum rules and the high estimate of Buch- 
mann [a b ased on a pion cloud model and the GP formal- 
ism I, 111 [3. As we have previously emphasized, the 
small result for 0^+ obtained from Sp 1 shows the im- 
portance of using the lattice data to constrain the model; 
without this constraint the uncertainty in our prediction 
of 0^+ would be much larger. 

Using the "minimal electromagnetic current" defined 
in the historical literature, these results imply an oblate 
form for both charge and magnetic distributions of the 
A^". However, the electromagnetic coupling recently de- 



4 



scribed by Alexandrou et al. [l7|, predicts that even a 
point-like A will have "natural" moments of Qa+ = 
-0.077 (Gfa (0) = -3) and Oa+ = -0.0021 (6^^3(0) = 
— 1) leading to a different interpretation of our results. 
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